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We report a direct observation of radiation pressure, 
exerted on cold rubidium atoms while bouncing on an 
evanescent- wave atom mirror. We analyze the radiation pres- 
sure by imaging the motion of the atoms after the bounce. 
The number of absorbed photons is measured for laser de- 
tunings ranging from 190 MHz to 1.4 GHz and for angles 
from 0.9 mrad to 24 mrad above the critical angle of total 
internal reflection. Depending on these settings, we find ve- 
locity changes parallel with the mirror surface, ranging from 
1 to 18 cm/s. This corresponds to 2 to 31 photon recoils per 
atom. These results are independent of the evanescent- wave 
optical power. 
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I. INTRODUCTION 

An evanescent wave (EW) appears whenever an elec- 
tromagnetic wave undergoes total internal reflection at a 
dielectric interface. The EW is characterized by an elec- 
tric field amplitude that decays exponentially with the 
distance to the interface. The decay length is on the or- 
der of the (reduced) optical wavelength. Cook and Hill 
[0 proposed to use the EW as a mirror for slow neutral 
atoms, based on the "dipole force". EW mirrors have 
since become an important tool in atom optics . They 
have been demonstrated for atomic beams at grazing in- 
cidence H and for ultracold atoms at normal incidence 

Most experimental work so far, has been concentrated 
on the reflective properties, i.e. the change of the atomic 
motion perpendicular to the surface [||. This is domi- 
nated by the dipole force due to the strong gradient of the 
electric field amplitude. In the present paper we report 
on our measurement of the force parallel to the surface. 
It has been mentioned already in the original proposal 
[0 that there should be such a force. The propagating 
component of the wavevector leads to a spontaneous scat- 
tering force, "radiation pressure" Hf?J. We present here 
the first direct observation of radiation pressure exerted 
by evanescent waves on cold atoms. Previously, a force 
parallel to the surface has been observed for micrometer- 
sized dielectric spheres ||. 

In our experiment, we observe the trajectory of a cloud 
of cold rubidium atoms falling and bouncing on a hori- 
zontal EW mirror. The radiation pressure appears as 



a change in horizontal velocity during the bounce. We 
study the average number of scattered photons per atom 
as a function of the detuning and angle of incidence of 
the EW. The latter varies the "steepness" of the optical 
potential. 

Due to its short extension, on the order of the opti- 
cal wavelength, an EW mirror is a promising tool for 
efficient loading of low-dimensional optical atom traps 
in the vicinity of the dielectric surface ]To],[ll|. In these 
schemes, spontaneous optical transitions provide dissipa- 
tion and cause the mirror to be "inelastic" , such that the 
final atomic phase- space density increases. This might 
open a route towards quantum degenerate gases, which 
does not use evaporative cooling and could poten- 
tially yield "atom lasers" ]l3| which are open, driven sys- 
tems out of thermal equilibrium, similar to optical lasers. 
It is this application of EW mirrors which drives our in- 
terest in experimental control of the photon scattering of 
bouncing atoms. 

This article is structured as follows. In Sec. II we sum- 
marize the properties of the EW potential, including the 
photon scattering of a bouncing atom. In Sec. Ill we 
describe our experimental setup and our observation of 
radiation pressure. We investigate its dependence on the 
angle of incidence and laser detuning. Finally, we discuss 
several systematic errors. 

II. THE EVANESCENT- WAVE ATOMIC MIRROR 

A. Optical dipole potential 

Our evanescent wave is created by total internal re- 
flection at a glass surface in vacuum Q . We choose the 
z-direction as the surface normal and the xz plane as the 
plane of incidence (see Fig. |l|). The EW can be writ- 
ten as E (r) oc exp(ik • r), where k = (k x ,0,iK). The 
wavevector has a propagating component along the sur- 
face, k x = /consin#, where ko = 2tt/\o is the vacuum 
wave number, n is the refractive index and the angle 
of incidence. Note, that k x > &o, since 6 is larger than 
the critical angle C = arcsinn -1 . The wavevector com- 
ponent perpendicular to the surface is imaginary, with 

k = kov n 2 sin 2 — 1. 
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FIG. 1. Evanescent- wave atom mirror, (a) Configuration 
in the rubidium vapor cell: magneto-optical trap (MOT), 
right-angle prism with refractive index n (6.6 mm below the 
MOT, gravity 1 1 2), evanescent- wave beam (EW), camera fac- 
ing from the side (CCD, in y direction), resonant fluorescence 
probe beam from above (P). (b) Confocal relay telescope for 
adjusting the angle of incidence 0. The lenses Li s 2 have equal 
focal length, / = 75 mm. A translation of L\ by a distance 
Aa changes the angle of incidence by AO = Aa/fn. The 
position of the EW spot remains fixed. M is a mirror. 

The optical dipole potential for a two-level atom at 
a distance z above the prism can then be written as 
Mdip(z) = Uq exp(— 2kz). In the limit of large laser de- 
tuning, \S\ ^> T, and low saturation, so < 1, the maxi- 
mum potential at the prism surface is Uq = h5so/2, with 
a saturation parameter so — (T/25) 2 TI/Iq [ fT5| . Here, 
Io = 1.65 mW/cm 2 is the saturation intensity for rubid- 
ium and T = 2tt x 6.0 MHz is the natural linewidth. The 
intensity of the incoming beam in the glass substrate is 
given as /. It is enhanced by a factor T, that ranges be- 
tween 5.4 and 6.0 for our TM polarized EW @. The 
detuning of the laser frequency uol with respect to the 
atomic transition frequency c^o is defined as 5 = uol — ooq. 
Thus a "blue" detuning, 5 > 0, yields an exponential 
potential barrier for incoming atoms. Given an incident 
atom with momentum p^, a classical turning point of the 
motion exists if the barrier height exceeds the kinetic en- 
ergy pf/2M of the atom. 

For a purely optical potential, the barrier height is 
given by Uq. In reality, the potential is also influenced by 
gravity and the attractive van der Waals potential, 



The gravitational potential U grav (z) oc z can be neglected 
on the length scale of the EW decay length £ = 1/ n. The 
van der Waals potential U v dw(z) cx (koz)~ 3 significantly 
lowers the maximum potential close to the prism and 
thus decreases the effective mirror surface on which atoms 
still can bounce. This effect has been demonstrated by 
Landragin et al Jlj. 



B. Photon scattering by bouncing atoms 

The photon scattering rate of a two-level atom in 
steady state at low saturation can be written as V = 
sT/2 = (T/hS)Udi P [15 1- An atom bouncing on an EW 
mirror sees a time dependent saturation parameter s(t). 
Assuming that the excited state population follows adi- 
abatically, we can integrate the scattering rate along an 
atom's trajectory to get the number of scattered photons, 
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For a purely optical potential, U cx exp(— 2k,z), this leads 
to an analytical solution: 
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The "steepness" of the optical potential is determined 
by n. The steeper the potential, the shorter the time an 
atom spends in the light field, and the lower N scat . Note 
that As Ca £ is independent of as a consequence of the 
exp(— 2kz) shape of the potential. This is in fact the rea- 
son why a two-level description of the atoms is appropri- 
ate. Strictly speaking, for a realistic atom the potential 
strength depends on the magnetic sublevel mp through 
the Clebsch-Gordan coefficients. For example, our EW 
is TM-polarized which, close to the critical angle, is ap- 
proximately linear. Effectively, for every mp sublevel, Uq 
is then multiplied by the square of a Clebsch-Gordan co- 
efficient, ranging from 1/3 to 3/5 (for aF = 2^F' = 3 
transition). Remarkably though, the value of A^ sca ^ re- 
mains unaffected and the same for all atoms. 

One expects that an absorbed photon gives a recoil 
momentum p rec = fik x to the atom, directed along the 
propagating component of the EW. Experimentally we 
observe this effect by the altered horizontal velocity of 
atom clouds after the bounce. The spontaneous emis- 
sion of photons leads to heating of the cloud and thus to 
thermal expansion |l8| , |l9| ]. 

In principle, A^ sca ^ is changed if other than optical 
forces are present. For example, the van der Waals at- 
traction tends to "soften" the potential and thus to in- 
crease N scat . We investigated this numerically and found 
it to be negligible in our present experiment. 
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III. THE OBSERVATION OF RADIATION 
PRESSURE 

A. Experimental setup 

Our experiment is performed in a rubidium vapor cell. 
We trap about 10 T atoms of 87 Rb in a magneto-optical 
trap (MOT) and subsequently cool them in optical mo- 
lasses to 10 fiK. The MOT is centered 6.6 mm above the 
horizontal surface of a right-angle BK7 prism (n = 1.51, 
6 C = 41.4° ||]), as shown in Fig. 0a. 

The EW beam emerges from a single-mode optical 
fiber, is collimated and directed to the prism through 
a relay telescope (see Fig.jjjb). The angle of incidence 
is controlled by the vertical displacement Aa of the first 
telescope lens, L\. This lens directs the beam, whereas 
the second lens, L2, images it to a fixed spot at the prism 
surface. A displacement Aa leads to a variation in 0, 
given by AO = Aa/nf. The focal length of both lenses is 
/ = 75 mm. The beam has a minimum waist of 335 /im 
at the surface (1/e 2 intensity radius). We checked the 
beam collimation and found it almost diffraction limited 
with a divergence half-angle of less than 1 mrad. We use 
TM polarization for the EW because it yields a stronger 
dipole potential than a TE polarized beam of the same 
power. Close to the critical ang le, the ratio of potential 
heights is approximately n 2 |1(| • 

For the EW, an injection- locked single mode laser 
diode (Hitachi HL7851G98) provides us with up to 
28 mW of optical power behind the fiber. It is seeded by 
an external grating stabilized diode laser, locked to the 
87 Rb hyperfine transition 5S 1/2 (F = 2) 5P 3/2 (F' = 3) 
of the D 2 line (780 nm), with a natural linewidth T = 
2tt x 6.0 MHz. The detuning 5 with respect to this transi- 
tion defines the optical potential Udi P for atoms that are 
released from the MOT in the F — 2 ground state. We 
adjust detunings up to 2tt x 500 MHz by frequency shift- 
ing the seed beam from resonance, using an acousto-optic 
modulator. For larger detunings we unlock the seed laser 
and set its frequency manually, according to the reading 
of an optical spectrum analyzer with 1 GHz free spectral 
range. 

Atoms that have bounced on the EW mirror are de- 
tected by induced fluorescence from a pulsed probe beam, 
resonant with the F = 2 — » F' = 3 transition. The probe 
beam travels in the vertical downward direction and has 
a diameter of 10 mm. The fluorescence is recorded from 
the side, in the y direction, by a digital frame-transfer 
CCD camera (Princeton Instruments). The integration 
time is chosen between 0.1 ms and 1 ms and is matched 
to the duration of the probe pulse. Each camera image 
consists of 400 x 400 pixels, that were hardware-binned 
on the CCD chip in groups of four pixels. The field of 
view is 10.2 x 10.2 mm 2 with a spatial resolution of 51 /im 
per pixel. 

A typical timing sequence of the experiment is as fol- 
lows. The MOT is loaded from the background vapor 



during 1 s. After 4 ms of polarization gradient cool- 
ing in optical molasses the atoms are released in the 
F = 2 ground state by closing a shutter in the cool- 
ing laser beams. The image capture is triggered with 
a variable time delay after releasing the atoms. During 
the entire sequence, the EW laser is permanently on. In 
addition, a permanent repumping beam, tuned to the 
5S 1/2 (F = 1) 5Pi /2 (F / = 2) transition of the D ± 
line (795 nm), counteracts optical pumping to the F = 1 
ground state by the probe. We observed no significant 
influence on the performance of the EW mirror by the 
repumping light. 

We measure the trajectories of bouncing atoms by tak- 
ing a series of images with incremental time delays. A 
typical series with increments of 10 ms between the im- 
ages is shown in Fig. |[ Our detection destroys the atom 
cloud, so a new sample was prepared for each image. The 
exposure time was 0.5 ms. Each image has been averaged 
over 10 shots. The fourth image shows the cloud just 
before the average bouncing time, tb = 36.7 ms, corre- 
sponding to the fall height of 6.6 mm. In later frames we 
see the atom cloud bouncing up from the surface. Close 
to the prism, the fast vertical motion causes blurring of 
the image. Another cause of vertical blur is motion due 
to radiation pressure by the probe pulse. The horizon- 
tal motion of the clouds was not affected by the probe. 
We checked this by comparing with images taken with 
considerably shorter probe pulses of 0.1 ms duration. 
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FIG. 2. Fluorescence images of a bouncing atom cloud. 
The first image was taken 5 ms after releasing the atoms from 
the MOT. The configuration of prism and evanescent wave is 
illustrated by the schematic (Field of view: 10.2 x 10.2 mm 2 ). 



B. Results 

Radiation pressure in the evanescent wave was ob- 
served by analyzing the horizontal motion of the clouds. 
From the camera images we determine the center of mass 
position of the clouds to about ±1 pixel accuracy. In 
Fig. [3], we plot the horizontal position vs. time elapsed 
since release from the MOT. We find that the horizon- 
tal motion is uniform before and after the bounce. The 
horizontal velocity changes suddenly during the bounce 
as a consequence of scattering EW photons. The change 
in velocity is obtained from a linear fit. 
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FIG. 3. Horizontal motion of bouncing atom clouds. The 
center of mass position is plotted vs. time elapsed since 
release. The bounces occur at 36.7 ms (vertical dotted 
line). (a) The EW decay length is varied as £ (0) — 
1.87, 1.03, 0.79, 0.67, 0.59, 0.53 fim (from large to small 
change in velocity). The detuning is 44 T and the optical 
power is 19 mW. (b) Comparison of two values of EW optical 
power, 19 mW (solid points) and 10.5 mW (open points). The 
detuning is 31 T and the EW decay lengths are £ = 1.87 /im 
(large velocity change) and 0.67 /im (small change). Solid 
lines indicate linear fits. 
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This corresponds with an EW decay length £(0) = 2.8 /im 
and 0.67 /im, respectively. We find that the number of 
scattered photons is inversely proportional to 5, as ex- 
pected. The predictions based on Eq. (||) are indicated 
in the figure (solid lines). 

In Fig.||b, the detuning was kept fixed at 44 T and 
the angle of incidence was varied between 0.9 mrad and 
24.0 mrad above the critical angle 9 C . This leads to a 
variation of £ (0) from 2.8 /im to 0.53 /im. Here also, we 
find a linear dependence on £ (0) . We see clearly, that a 
steep optical potential, i.e. a small decay length, causes 
less radiation pressure than a shallow potential. 

The observed radiation pressure ranges from 2 to 31 
photon recoils per atom. Note, that we separate this 
subtle effect from the faster vertical motion, in which 
atoms enter the optical potential with a momentum of 

Pi ~ 61 Prec 

In Fig.[|b, we compare trajectories for 19 ± 1 mW and 
10.5 ± 0.5 mW optical power in the EW. As expected 
from Eq. (||), there is no significant difference in horizon- 
tal motion. For a decay length of £ = 2.78 /im (0.67 /im) 
both power settings lead to essentially the same radia- 
tion pressure, that is 25 ± 3 (13 ± 2) scattered photons 
for 19 mW, and 23 ± 2 (11 ± 1) photons for 10.5 mW. 
The optical power only determines the effective mirror 
surface and thus the fraction of bouncing atoms. This is 
also visible in the horizontal width of bouncing clouds. 
The bouncing fraction scales with the intensity and the 
detuning as oc \n(I/S). We observe typical fractions of 
13 % for 5 = 44 T. For a given optical power, there is an 
upper limit for the detuning, above which no bounce can 
occur. For the data in Fig.[|a, the threshold is calculated 
as 6th = 6.5 GHz (8.1 GHz) for £ = 2.8 /im (0.67 /im). 
Another threshold condition is implied by the van der 
Waals interaction, and yields a lower limit for the EW 
decay length £. For Fig. [|b this lower limit is calculated 
as ^th = H6 nm, i.e. t h = (0 C + 0.59 rad). 



C. Systematic errors and discussion 



FIG. 4. Radiation pressure on bouncing atoms expressed 
as number of absorbed photons, J\f sca t- (a) Detuning 5 varied 
for £ = 2.8 /im (open points) and 0.67 /im (solid points), (b) 
EW decay length £ varied for 5 — 44 T. The laser power 
was 19 mW. The thin solid line is a linear fit through the first 
four data points. Theoretical predictions: two-level atom (see 
Eq.(^|), thick solid lines). Rubidium excited-state hyperfine 
structure and saturation taken into account (dashed lines). 

In Fig.||, we show how the radiation pressure depends 
on the laser detuning 5 and on the angle of incidence 0. 
The fitted horizontal velocity change has been expressed 
in units of the EW photon recoil, p rec = hkonsinO, with 
hko/M 5.88 mm/s. 

In Fig.||a, the detuning is varied from 188-1400 MHz 
(31 — 233 T). Two sets of data are shown, taken for two 
different angles, = 6 C + 0.9 mrad and C + 15.2 mrad. 



According to Eq. (||), the radiation pressure should be 
inversely proportional to both 5 and k(Q). As shown in 
Fig. [|, we find deviations from this expectation in our ex- 
periment, particularly in the k dependence. A linear fit 
to the data for £ < 1 /im, extrapolates to an offset of ap- 
proximately 3 photon recoils in the limit £ — » (thin solid 
line in Fig.^b). The vertical error bars on the data in- 
clude statistical errors in the velocity determination from 
the cloud trajectories as well as systematic errors. We 
discuss several possible systematic errors, namely (i) the 
geometric alignment, (ii) the EW beam angle calibration 
and collimation, (iii) diffusely scattered light, (iv) the 
van der Waals atom-surface interaction, (v) excited state 
contributions to the optical potential, (vi) and saturation 
effects. 

(i) Geometrical misalignments give rise to systematic 
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errors in the radiation pressure measurements. For ex- 
ample, a tilt of the prism causes a horizontal velocity 
change even for specularly reflected atoms. We checked 
the prism alignment and found it tilted 12 ±5 mrad from 
horizontal. This corresponds to an offset of 1.5 ± 0.6 re- 
coils on N scat . In addition, the atoms are "launched" 
from the MOT with a small initial horizontal velocity 
which we found to correspond to less than ±0.4 recoils 
for all our data. From Fig.||, we see that the extrapolated 
trajectories at the bouncing time do not start from the 
horizontal position before the bounce. We attribute this 
to a horizontal misalignment of the MOT with respect to 
the EW spot. Obviously, there is a small displacement of 
the EW spot at the prism surface, when adjusting 6 by 
means of the lens L\ (see Fig.|l|b). Since the finite sized 
EW mirror reflects only part of the thermally expand- 
ing atom cloud, such a displacement selects a nonzero 
horizontal velocity for bouncing atoms. We corrected for 
those alignment effects in the radiation pressure data of 
Fig.[|. For small radiation pressure values, the system- 
atic error due to alignment is the dominant contribution 
in the vertical error bar. 

(ii) The uncertainty in the EW angle with respect to 
the critical angle is expressed by the horizontal error 
bars. We determined — C within ±0.2 mrad by mon- 
itoring the optical power transmitted through the prism 
surface, while tuning the angle from below to above 
6 C . Close to the critical angle, the decay length £(0) di- 
verges, and thus the error bar on £ becomes very large. 
Also the diffraction-limited divergence of the EW beam 
may become significant. It causes part of the optical 
power to propagate into the vacuum. In addition, the 
optical potential is governed by a whole distribution of 
decay lengths. Thus the model of a simple exponential 
optical potential oc exp {—2nz) might not be valid and 
cause the disagreement of our data with the prediction 
by Eq. (H). For larger angles, i.e. £(#) < 1 /im, the effect 
of the beam divergence is negligible. This we could verify 
by numerical analysis. 

(hi) Light from the EW can diffusely scatter and prop- 
agate into the vacuum due to roughness of the prism 
surface. We presume this is the reason for the extrapo- 
lated offset of 3 photon recoils in the radiation pressure 
(Fig. Qd). A preferential light scattering in the direction 
of the EW propagating component can be explained, if 
the power spectrum of the surface roughness is narrow 
compared to 1/ A (l9) . The effect of surface roughness on 
bouncing atoms has previously been observed as a 
broadening of atom clouds by the roughness of the dipole 
potential. In our case, we observe a change in center 
of mass motion of the clouds due to an increase in the 
spontaneous scattering force. Such a contribution to the 
radiation pressure due to surface roughness vanishes in 
the limit of large detuning S. Thus, we find no significant 
offset in Fig.[|a. Scattered light might also be the rea- 
son for the small difference in radiation pressure for the 
two distinct EW power settings, shown in Fig.||b. Lower 
optical power implies slightly less radiation pressure. 



(iv) As stated above, the van der Waals interaction 
"softens" the mirror potential. This makes bouncing 
atoms move longer in the light field, thus enhancing pho- 
ton scattering. We investigated this numerically by inte- 
grating the scattering rate along the atoms' path, includ- 
ing the van der Waals contribution to the mirror poten- 
tial. Even for our shortest decay parameter of 0.53 /im, 
the number of scattered photons would increase only 
about 0.8% compared with Eq. ([|), which we cannot re- 
solve experimentally pl| . 

(v) In a two-level atom model, the scattering rate can 
be expressed in the dipole potential as r' = (Y /ft5)Udi p - 
This is no longer true if we take into account the excited 
state manifold F' = {0, 1,2,3} of 87 Rb. Beside F' = 3, 
also F' = 2 contributes significantly to the mirror poten- 
tial, whereas it does not much affect the scattering rate. 
With an EW detuning of 44 T, this results in a number 
of scattered photons, N scat , typically 9% lower than ex- 
pected for a two- level atom, Eq. (||). Here we averaged 
over contributions from distinct magnetic sublevels. 

(vi) In order to investigate the influence of saturation 
on the number of scattered photons, N scatl we solved 
the optical Bloch equations numerically. A bouncing 
atom encounters the EW as a fast varying light pulse 
I(t) oc sech 2 (ftpit I 'M) , with a typical duration between 
3 and 10 fis. By integrating the time-dependent scatter- 
ing rate for an atom bouncing with an EW detuning of 
44 T, we find approximately 7% less scattered photons 
compared with the unsaturated expression of Eq. (||). 
Note, that the bounces occur sufficiently slow to pre- 
serve adiabaticity. In Fig.[|, we show predicted curves, 
corrected for hyperfine structure and saturation (dashed 
solid lines). 

IV. CONCLUSIONS 

We have directly observed radiation pressure that 
is exerted on rubidium atoms while bouncing on an 
evanescent- wave atom mirror. We did so by analyzing 
the bouncing trajectories. The radiation pressure is di- 
rected parallel to the propagating component of the EW, 
i.e. parallel to the interface. We observe 2 to 31 photon 
recoils per atom per bounce. We find the radiation pres- 
sure to be independent of the optical power in the EW, 
as expected from the exponential character of the EW. 

The inverse proportionality to both the EW detuning 
and the angle of incidence is in reasonable agreement 
with a simple two-level atom calculation, using steady 
state expressions for the EW optical potential and the 
photon scattering rate. The agreement improved when 
also the excited state hyperfine structure and saturation 
effects were taken into account. The measured number 
of photon recoils as a function of decay length £ indicates 
an offset of approximately 3 recoils in the limit of a very 
steep EW potential. We assume this is due to light that is 
diffusely scattered due to roughness of the prism surface, 
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but retains a preferential forward direction parallel with 
the EW propagating component. 

With sufficient resolution, it should be possible to re- 
solve the discrete nature of the number of photon recoils 
and also their magnitude, Tik x > hko |J. Our technique 
could also be used to observe quantum electrodynamical 
effects for atoms in the vicinity of a surface, such as ra- 
diation pressure out of the direction of the propagating 
EW component pgfl. 
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